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Abstract. I show how transition systems can be applied to the nat- 
urally concurrent behaviour of excitable media. I consider structured 
excitable media, in which excitations are constrained to propagate only 
in defined narrow channels, and cannot propagate elsewhere. I define a 
type of transition system that can be used to describe the complete set 
of behaviours exhibited by simple structures. The composition rules that 
result from this definition can be used to automatically deduce the be- 
haviour of more complex structures composed from simpler structures. 
Several examples illustrate the method, and a software implementation 
is provided. 
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1 Introduction 

Unconventional computing seeks to explore and characterise a wider range of 
substrates and phenomena than those that have been used so far for practical 
computing, with the hope of understanding and applying the intrinsic informa- 
tion processing behaviour of each substrate. Unconventional computing tends to 
be substrate-driven rather than model-driven. Rather than searching for physi- 
cal systems that can be harnessed to make computing machinery according to 
a preconceived idea of which behaviour the machinery is required to exhibit, 
unconventional computing is inspired by the idea that every sufficiently complex 
system performs information processing of some kind, which we can harness if 
we can find the appropriate mathematical techniques for understanding or in- 
terpreting the system, and the appropriate experimental techniques for manipu- 
lating and observing the system. The existence of natural systems that manage 
to harness complex phenomena to process information in powerful ways that 
we do not yet fully understand is one of the motivations for this search, though 
unconventional computing does not confine itself to naturally occurring systems. 

One of the difficulties in efficiently exploiting previously unexploited phe- 
nomena in new or existing substrates is that many of the abstract models that 
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we use in computer science were invented and popularised precisely because they 
abstract away from the details of the physical behaviour of computing hardware. 
In the early history of computing machinery it was often pragmatic to regard 
concurrent, time-dependent behaviour as a nuisance to be abstracted away from 
so that it did not complicate and impede design. In synchronous circuit design 
the input-output delay present in all logic gates is taken account of by waiting 
for a sufficiently long time that all gates have had enough time to propagate any 
changes on their inputs to their outputs, and only then carrying out the next 
step in the computation. Asynchronous design techniques [1] can help designers 
to make more efficient designs by having circuits explicitly signal when their 
outputs become valid, so that there is no need to wait for the worst-case delay 
time, and no need to have all circuits in a system clocked at constant rates. 

Asynchronous circuit design involves understanding electronic digital circuit 
elements at a less abstract level than that which is used in synchronous design, 
this makes it possible to exploit types of behaviour that do not exist at higher 
levels of abstraction. This approach towards the substrate can be generalised 
[2] ■ Some of the questions that can be asked of any substrate considered to have 
complex information processing potential are: 

1. What are the intrinsic behaviours of a region of substrate and how can they 
be modelled? 

2. How do regions of substrate behave when connected together? 

3. Can the modes of behaviour of regions of substrate, or composed regions of 
substrate be mapped easily onto any existing abstractions used for informa- 
tion processing? 

4. Does the substrate suggest new models for information processing? 

5. Can the substrate be applied to a problem without having to have a detailed 
abstract model for understanding how it behaves as it does? 

For complex substrates answers to questions 3,4 and 5 can be sought through 
observation and experimentation without having quantitative answers for ques- 
tions 1 and 2. For example, the tendency of plasmodium of physarum poly- 
cephalum to form efficient networks connecting multiple food sources was jus- 
tified on evolutionary grounds and observed experimentally [3] before a math- 
ematical model of dynamic network formation was constructed [4], and before 
a low-level particle based model that can give rise to that behaviour was con- 
structed [5], Robots have been devised that have controllers which exploit a 
substrate without requiring a theory that accounts for the behaviour of the sub- 
strate [5T7] , 

For other substrates, the behaviour of small regions of substrate may be 
easy to explain and describe [5], but challenges arise with recognising that the 
substrate might be capable of information processing, with mapping it onto 
higher level abstractions, and with trying to deduce the behaviour of large regions 
of substrate from the behaviour of small regions. 

In this paper I consider excitable media substrates, which can be used to 
implement channels along which excitations can propagate. Examples of this 
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type of substrate are: planar Belousov Zhabotinsky (BZ) reaction media, top- 
pling dominoes, plasmodium of physarum polycephalum, and the propagation 
of flames in a forest fire. Excitable media are characterised by the property that 
waves of disturbance can propagate through a medium without diminishing. The 
physical phenomenon that gives rise to this behaviour differs from one medium 
to another. The toppling domino and the forest fire substrates are examples of 
single-use media — a disturbance can propagate through a region once only. 
For the sake of simplicity I consider only structured excitable media, where ex- 
citations are confined to predefined regions, and do not consider unstructured 
free-space systems, where excitations propagate and interact in a homogeneous 
medium [5]. 

One way of using excitable media for information processing is to exploit the 
property that some media have whereby two travelling excitations close to each 
other can exert an influence over each other. The consequences of that influence 
can be detected by exploiting another property: relative to their direction of 
movement, excitations may spread more favourably in some directions than in 
others. Flames propagating in a forest fire tend to spread out in all directions 
- a flame front passing by a region of unburnt forest will not leave that part 
unburned. In contrast, a toppling domino will only topple in a direction per- 
pendicular to its base. A toppling domino will not cause its neighbour to topple 
unless that neighbour happens to stand in the way of the toppling domino. In 
between these two extremes we find other types of propagation behaviour which 
can favour one locus over another. Planar BZ-reaction media can be made to 
exhibit directionally dependent propagation in a number of ways jlOll lj . The 
preference of plasmodium of physarum polycephalum for one direction over an- 
other can be influenced by many different factors |12) . In certain circumstances 
the behaviour of propagating plasmodium has striking similarities to the be- 
haviour of propagating waves in BZ reaction media [13114] and can similarly be 
made to exhibit directionally dependent propagation (Chapter 7 of |12jV 

A recent survey of some of the ways in which structured BZ-reaction media 
have been used for information processing is given by Gorecki |llj . Excitable 
media have been simulated at a low level of abstraction using partial differential 
equation models, and these low-level models have been used (in conjunction with 
experiments) to gain qualitative insight into the behaviour of the media and to 
discover relatively high-level structures such as Boolean logic gates which can be 
put together to make complex circuits. I attempt to fill a gap between these two 
levels of abstraction by quantitatively describing the time-dependent behaviour 
of regions of structured excitable media in response to excitations. Describing 
structures at this level of abstraction retains timing information, which can be 
exploited to make more efficient circuits, or to make systems with behaviour 
that cannot be succinctly described in terms of combinational or sequential logic 
circuits. 

If a pattern of excitation is applied to an excitable medium it will evolve 
over time in a generally deterministic way. If another excitation is applied to the 
medium as it is evolving, the way in which it evolves may change, and it might 
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unfold into a different sequence of patterns. In general, the way that the system 
evolves over time is dependent on the timing and location of the excitations that 
it receives. Within concurrency theory the term reactive system is often used 
to describe a system whose behaviour depends on the pattern and timing of 
events that influence the system [TS]. I would like to avoid potential confusion 
by making it clear at this point that the term reactive system has nothing to do 
with the fact that some excitable media make use of chemical reactions. 
In summary, the two problems that I address in this paper are: 

1. How can regions of excitable media be described as reactive systems'! 

2. How can the behaviour of complex structures made from excitable media be 
deduced from the behaviour of their constituent structures? 

This work has similar aims to some of the uses of concurrency techniques in 
asynchronous circuit design |16) . One difference is that in asynchronous design 
these techniques are used for the specification and verification of relatively large 
systems whereas here the emphasis is on exploring and uncovering the dynamics 
of relatively small structures. Another difference is that in asynchronous systems 
few assumptions are made about timing, but here all signals are assumed to 
propagate at a known speed. 

The purpose of this work is not simply to simulate the behaviour of structures 
made from excitable media at an abstract level — this can be done perfectly 
well using cellular automaton models of excitable media. Rather, this work is 
concerned with describing the complete set of behaviours that a structure can 
exhibit, so that these behaviours can be understood and reasoned about without 
having to run exhaustive simulations of a system. 

2 Fork structures in excitable media 

Figure [T] shows a domino fork. This is a useful structure to use for the expo- 
sition which follows because it is a simple, familiar system. It is arguably the 
simplest structure in an excitable medium that can be used as the sole basis for 
more complex information processing structures such as Boolean logic gates and 
circuits. 

The behaviour of the domino fork is summarised in the following paragraphs. 
Here the word 'event' is used to mean the toppling of a domino. Dominoes can 
topple in any one of two directions. An 'input event' to a system of dominoes 
is a domino outside the system toppling towards the system. An 'output event' 
is a domino inside the system toppling away from and out of the system. The 
word 'signal' is used to mean a travelling wave of excitation. 

The domino fork has three regions on its boundary in which events can occur. 
In each region there are two possible events that can occur: an input event, and 
an output event. We use a different label for each event. Clearly in a system 
of dominoes once one event has happened, the other event cannot happen. In 
reusable excitable media, which can propagate excitations in the same region 
again and again, this is not the case. The definitions of 'input event' and 'output 



Transition systems for excitable media 



5 



a 




Fig. 1. A domino fork. 



Fig. 2. A domino line. 



event' given above permit an input event and output event in the same region to 
occur at the same time: it is possible for a domino outside the system to topple 
towards the system at the same instant that a domino inide the system topples 
away from the system. The two dominoes will meet each other in the middle of 
toppling. 

If event c occurs then events d and / will occur a short time later, unless 
events a or e occur in the meantime. 

If event a or event e occur then event b will occur a short time later, unless 
event c occurs in the meantime. 

If we ignore one of the two branches of a fork, we end up with a simple line 
of dominoes, shown in Figure [2j where input event a will cause output event b a 
short time later, unless event c occurs in the meantime (with identical behaviour 
in the other direction). 

O'Keefe [T7] showed that, subject to certain timing constraints and assuming 
that bridges can be built to cross one line of dominoes over another, the domino 
fork can be used to build Boolean circuits. I showed that any Boolean circuit 
can be built, that bridges are not required, and that an asynchronous scheme 
can be used in which there are no timing constraints [18] . albeit with a larger 
number of forks for a given circuit than in O'Keefe's scheme. 

A domino fork can be used as a direction-detecting gate. To help conceptu- 
alise this, imagine bending the arms of the fork into a T shape (Figure [3]), and 
regard the horizontal channel as a channel in which a signal may propagate in 
either direction, while the vertical channel is treated as an output. Now a signal 
travelling from right to left in the horizontal channel will cause an output to 
emerge from the vertical channel, but a signal travelling from left to right will 
not result in an output from the vertical channel. 

The behaviour of a domino line can be regarded from the perspective of 
collision-based computing |19j . Here, a domino line represents a channel into 
which objects — in this case toppling wavefronts — can be sent. If two objects 
are sent into the channel in opposite directions they will collide and annihilate 
with each other, but if only one object is sent in it will emerge from the other 
end of the line. Therefore a line effects an AND-NOT operation: if event a 
AND-NOT event c then event b (subject to timing constraints). Domino fork 
structures placed at either end of a line can be seen as transducers which control 
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Fig. 3. The domino fork as a direction detecting gate. 



the injection of objects into the line, and which catch objects emerging from the 
line. 

Both the AND-NOT behaviour of a domino line and the direction-detecting 
behaviour of a fork seem to be essential for being able to construct more complex 
information processing structures using forks connected by lines. 

In [IS], several simple structures made from domino forks and lines were 
described, before it was shown how they can be put together to make arbitrary 
Boolean circuits. Two of these simple structures are shown in Figures 2] and [5] 



B" A' 




Fig. 4. A one way line. Fig. 5. A single-line crossover. 

Figure |4] shows a one way line constructed using two forks. A signal entering 
the configuration from A will split at fork F into two signals that will collide with 
each other, preventing any signal from emerging at B. In the other direction, a 
signal entering the configuration from B will pass through fork G, through F 
and then emerge from A. 

Figure [5] shows a single-line crossover (so-called because it allows two single 
lines to cross over, whereas another structure in JT5] allows two pairs of lines to 
cross over). It has two inputs A and B and two outputs A' and B' which are 
crossed over topologically. A signal may arrive either at input A or at input B 
(but not both). A signal arriving at A will propagate to A', but not to B or B'. 
A signal arriving at B will propagate to B', but not to A or A'. 

The fork structure can also be implemented in other excitable media. Motoike 
and Yoshikawa give an example of a fork-like structure in simulated BZ-reaction 
media [TU] and use it for implementing an OR gate by making use of its be- 
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haviour for the cases when a or e or both receive excitations (with reference to 
Figure [I| . Motoike and Yoshikawa used unexcitable barriers to prevent back- 
ward propagation from one fork input to another, and as a consequence their 
fork effectively has no c input. An alternative way to implement a fork struc- 
ture is to use a light sensitive medium, in which the excitability of the medium 
can be controlled by the level of illumination [3D]. Within a certain range of 
illumination, the medium can be made to exhibit sub-excitable behaviour, in 
which excitations in channels cannot propagate around sharp corners. Figures [SJ 
and [7] show two snapshots from a simulation of a fork structure in this medium. 
Here the dark regions support sub-excitable waves, and the lighter surrounding 
regions inhibit excitation. This behaviour has also been observed experimentally 
in [H]. 




Fig. 6. Splitting at a fork junction in Fig. 7. Directional propagation at a 
sub-excitable BZ-reaction media. fork junction in sub-excitable BZ- 

reaction media. 



The behaviour of the fork structure in sub-excitable BZ-reaction media is 
very similar to the behaviour of the domino fork, and depends on the directional 
propagation of excitations in the sub-excitable regime. The most significant dif- 
ference between BZ-reaction media and domino media is that a BZ-reaction 
medium is reusable: once the medium has recovered from an excitation it can be 
used again. For some structures this can lead to a significantly different range 
of behaviours than those exhibited by a single use medium. This is discussed in 
section 14.21 

Plasmodium of physarum polycephalum has a mode of behaviour in which 
it moves along a regularly spaced line of food particles in a spreading and con- 
tracting mode [22] . If food particles are placed close together along a line then a 
Plasmodium will spread from a food particle until the front of the spreading area 
encounters another food particle. It will quickly slow and stop its spread, then 
migrate to the new food particle, leaving a tube connecting with the old parti- 
cle. It will then begin spreading again from the new particle, repeating the same 
behaviour. Figure [8] shows an image of a plasmodium exhibiting this behaviour. 
Figures [SJ and [TUJ show work in progress to coax physarum to exhibit fork-like 
behaviour, obtained from the spreading and contracting mode of propagation 
with a carefully designed food geometry. 
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Fig. 8. The spreading and contracting behaviour of Plasmodium of physarum poly- 
cephalum. 



Fork structures and line structures are used exclusively in the exposition 
that follows. The method that is introduced should work equally well for other 
structures with behaviour that can be described in terms of a reactive system. 

3 Transition systems 

I begin by defining a type of transition system that is capable of modelling 
the behaviour of single fork and line structures in terms of input and output 
events, and the timing relationships between them. I then give a set of rules for 
composing transition systems that will allow us to deduce the behaviour of larger 
and more complex structures from the behaviour of fork and line structures. 
When the total range of behaviours that a complex structure can exhibit is 
difficult to comprehend, I show how to select portions of its behaviour which can 
be understood in isolation and which can help to understand the total behaviour. 

The transition systems defined here make use of the description of transition 
systems given in [23J . The meaning given to transitions is similar to that used in 
SCCS [23] in so far as every transition represents the passage of a unit of time, 
in addition to the occurrence of zero or more events, but instead of using event 
sets as labels, Boolean expressions over event symbols are used. A common clock 
is assumed among all transition systems, so transitions in any pair of transition 
systems occur together at the same time. 




Fig. 9. Physarum exhibiting a splitting be- 
haviour at a fork junction. 



Fig. 10. Physarum exhibiting direction de- 
pendent propagation behaviour, induced by 
the food geometry. 
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Notation and concepts will be introduced by giving examples along with 
more formal descriptions. Whenever the term 'transition system' is used from 
here onwards, it refers to the particular type of transition system developed here, 
rather than to transition systems in general. 

Definition 1. A transition system is a tuple (S,i, L,Tran, E, I) where S is a 
set of states with initial state i, L is a set of labels, Tran C S x L x S is the 
transition relation, E is a set of event symbols, and I C E is a set of input event 
symbols. Symbols in E\I are output event symbols. 

If we are referring to two or more different transitions systems then in order 
to distinguish one from the other we use a dot notation to refer to the elements 
of a tuple: If P is a transition system then P.S is the set of states of P, P.i is 
the initial state of P, etc. 

Event symbols represent the occurrence of an event in the physical system 
that the transition system represents. For the sake of brevity we will use the 
word 'events' to refer to both event symbols and the events that the symbols 
represent. 

The set /, which we will call the input interface, is a set of events that are 
generated outside the transition system and which the system reacts to. The set 
E\I, which we will call the output interface, is a set of events that are generated 
by the system. 

A transition system represents a structure made from excitable media. Each 
state corresponds to a particular pattern of excitation, and transitions from one 
state to another correspond to the evolution over time of patterns of excita- 
tion, possibly subject to the effects of input events. The input interface of the 
transition system corresponds to locations where excitations can be fed into the 
structure, and the output interface corresponds to locations where excitations 
can emerge from the structure for observation or to feed into the inputs of other 
structures. 

Labels are Boolean expressions over events, where the symbol a in an ex- 
pression corresponds to the occurrence of event a, and the expression corre- 
sponds to the non-occurrence of event a. We will use Bool (A) to refer to the set 
of Boolean expressions over events from the set A. In our Boolean expressions 
we use • for conjunction, + for disjunction and -i for negation. Two labels are 
equivalent if they represent the same Boolean function, regardless of how that 
function is expressed. A label represents a collection of sets of events. We choose 
to use Boolean expressions (rather than event sets as in SCCS) for labels be- 
cause by doing so we often reduce the number of transitions that are needed in a 
transition system, because we can specify many sets of events in a single transi- 
tion, rather than having a separate transition for each set. The synchronisation 
operator for composing two transition system remains simple — it is essentially 
the conjunction of two Boolean expressions. 

Definition 2. We say that a label expression I with event variables aj matches 
a set of input events A C I if given the truth assignment aj = True if aj e 
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A, False if aj £ I \A, the resulting expression is not equal to False, i.e. it is 
either equal to True, or it has some remaining variables not in I. 

We must place some additional constraints on transition systems to bring 
their behaviour into closer agreement with the physical systems we are interested 
in. We define these constraints here, and later on when we compose two transition 
systems we will need to show that these constraints still hold in the transition 
system that is produced as the result of the composition. 

A transition corresponds to the passage of a unit of time in addition to 
the occurrence of events. We can think of the occurrence of events as being 
interleaved with the passage of units of time. Events occur, then a unit of time 
elapses, then more events occur, and so on. When a particular set of input events 
occurs, the transition followed is that which matches the events that occur. A 
transition system must be deterministic: no set of input events should match 
more than one transition from a given state. Additionally, we will adopt the 
convention that there must always be a transition that matches any given set of 
input events. Even if a system reaches a state from which it can never leave, we 
require that it has a transition back to the same state, matching any set of input 
events. The reason for this convention is that it simplifies the rule for composing 
two systems. 

Constraint 1 For each (s,A) 6 S x V{I) there must be exactly one (s,l,s') e 
Tran such that I matches A. 

A label must be the conjunction of two subexpressions, one being any Boolean 
expression of input events, the other being a conjunction of output events or their 
negation, in which every output event appears. This ensures that every transition 
effectively specifies which output events occur during that transition — they are 
the events that are not negated in the output event subexpression. 

We first define a set of expressions that are conjunctions of variables or their 
negation: 

Definition 3. Let Conj(A) = {b\ ■ b 2 ■ ... • b n \ (bi,b 2 , • ••,&«) £ n{ a ji"i}} f or 
aj e A 

Constraint 2 L = {x ■ y \ x e Bool(I) and y e Conj(E \I)} 

One of the advantages of using Boolean expressions representing collections 
of event sets for labels is that if a particular input event has absolutely no effect 
on the behaviour of a transition system in a given state — that is, the destination 
state is the same regardless of whether that input event occurs or not — then we 
can omit that event from any label in transitions from that state. A transition 
system in a state where some inputs have no effect could correspond to a physical 
system in a state where an input feeds into a region of excitable medium which 
has not yet recovered from a previous event. More generally it corresponds to a 
physical system evolving along a trajectory on which a particular input makes 
no difference to that trajectory. 
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The rule for determining which transition a transition system will undergo 
in response to a set of input events, and which output events it will produce is 
given by the following definition: 

Definition 4. If a system is in state s and the set of input events that occurs 
is A then the next state s nex t of the system is given by the element (s, I, s nex t) G 
Tran for which I matches A. The set of output events that the system will produce 
are those which appear non-negated in the output event subexpression of I. 

Definition 5. A state s is reachable if and only if it is the initial state i, or 
there is a transition (s', I, s) G Tran where I is satisfiable and s' is reachable. 

Constraint 3 All states in a transition system must be reachable. 

Finally, the next constraint specifies that no input can have an instantaneous 
effect on any output. There must be a delay of at least one time unit between 
the occurrence of a set of events and a change in the output behaviour of the 
system. This is realised by making sure that all labels for transitions from a 
given state have the same output event subexpression. 

Constraint 4 For each s G S there is a single expression y s G Conj(E \ I) 
such that for every (s, Z, s') G Tran, I is of the form x ■ y s (where x G Bool(I) is 
different for each transition). 

It is reasonable to insist that any transition system that corresponds to a 
physical excitable media system satisfies constraint |4l If we were to permit a 
transition system to flout constraint 2J then we could end up in a contradictory 
situation where two systems that shared both inputs and outputs could specify, 
say, that events a and b must always occur simultaneously, and also that b can 
only occur when a does not occur. 

These definitions and constraints permit us to specify transition systems that 
represent all of the possible behaviours of a fork and a line. Before doing so we 
will consider some simpler examples and then define rules for composing two 
transition systems. These simple examples do not correspond to any structures 
in excitable media. They serve simply to illustrate the concepts that have been 
described so far and introduce in an incremental way the notation and conven- 
tions used for representing transition systems graphically. 

Figure [11] is an example of a simple transition system. The input interface 
for this system is a and the output interface is b. This is displayed as a : b in 
the box in the bottom left of Figure QTJ The initial state is represented using 
one circle inside another and is also named as state 0. Other states are given 
numeric names. Transitions are represented as arrows from one state to another, 
with the label for the transition written above the arrow. 

This example introduces two conventions for representing transition systems 
that make them easier to read. The first convention is the way that label expres- 
sions are written. Labels could be written simply as Boolean expressions, but this 
representation cannot easily be read at a glance, so instead we take account of 
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a:b 



Fig. 11. An example of a transition system. 



constraint [5] regarding the allowable sub-expression for output events, and write 
the sub-expression corresponding to the input event variables first, followed by 
a colon, followed by a set of those output event variables that are not negated 
in the label expression. For example, if a and b are input events, and c and d are 
output events then the label expression (a + e)cd is written (a + e) : d. The colon 
is omitted if the output event subexpression is a conjunction of negated output 
event variables, or if the label expression matches any set of input events. If the 
label expression is simply a conjunction of negated output events then we write 
the label as 1. 

The second convention concerns the graphical representation of a consecu- 
tive sequence of transitions from one state to another that is not influenced by 
any input events and which does not produce any output events. Since such a 
sequence of states and transitions corresponds to the passage of a length of time 
equal to the length of the sequence, we write it using a single arrow with a num- 
ber above the arrow representing the length of the sequence. In Figure 111! the 
arrow with the number 5 written above it between states 1 and 2 has the same 
meaning as if we drew four extra states and transitions between state 1 and 2, 
and labelled each of the five transition between states 1 and 2 with the label 1. 

The behaviour of the transition system in Figure [11] can be summarised as 
follows. If the system is in its initial state and event a does not occur then it will 
remain in the initial state no matter how many time units elapse. If the system 
is in its initial state and input event a occurs then the system will be in state 
1 after 1 unit of time has elapsed. After a further 5 units of time have elapsed 
the system will enter state 2. It will then generate output event b and will be in 
state 3 after 1 unit of time. Once the system is in state 3 it will remain in state 
3 regardless of how much time elapses. 

Figure [T2] shows a slightly more complex example. The input interface of this 
system is a, & and the output interface is c,d,e. If the system is in its initial state 
and input event a occurs, or if both a and b occur, it will enter state 1 after 1 
unit of time. It will then generate event e and enter state 4 after a further unit of 
time. If the system is in its initial state and input event b occurs by itself, it will 
enter state 2 after 1 unit of time. It will then generate event d and enter state 4. 
If the system is in its initial state and neither input event occurs then one unit 
of time will elapse and the system will enter state 3. It will then output event 
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3.1 Composing transition systems 

In general, the behaviour of two transition systems P and Q operating inde- 
pendently and concurrently can be described by a single transition system by 
considering each state in P in combination with each state in Q. Figure [TBI shows 
two simple transition systems that have no events in common. Although these 
two transition systems operate independently, both are in their initial state at 
the starting time and they operate under a common clock. Because both of the 
transition systems in Figure [13] have transitions returning to every state, each 
of the systems can remain in the initial state or the final state for any length of 
time while the other system makes a state transition. Figure [T4l shows the result 
of composing these transition systems. 

The notation P \\ Q will be used for the composition of P and Q. When P 
and Q are composed the interface of P may have events in common with the 
interface of Q, so we must come up with a set of rules for combining transitions to 
make sure that the occurrence or non-occurrence of events in P is synchronised 
with the occurrence or non-occurrence of events in Q. 

For two transition systems to be composable, they must both satisfy con- 
straints [lj [2] and [3j and any event that they have in common must be an input 
in one system and an output in the other. 

Regarding constraint [H we will allow that either P or Q does not satisfy this 
constraint. The reason for this is made clear later on in section 14.21 However, 
whichever of P or Q does not satisify constraint 2] must have no inputs coming 
from the other system, and all of its outputs must feed into the other system. 

Definition 6. Two transition systems P and Q are composable if and only if 
the following conditions hold: 

A Both P and Q satisfy constraints^ [H andfJJ 
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Fig. 14. The composition of the two independent systems from Figure [131 
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B If P does not satisfy constraint ^ then Q must satisfy constraint [7] and 

P.E n Q.E = P.E \ P.I. 
C If Q does not satisfy constraint [^} then P must satisfy constraint and 

P.E n Q.E = Q.E\Q.I. 
Dae P.E n Q.E (a G P.I and a £ Q.I) or (a G Q.I and a £ P.I). 

To compose two transition systems P and Q to produce a transition system 
T = P || Q we first form the product T", where labels in transitions are combined 
by conjunction: 

Definition 7. The product T' = P x Q is defined as: 
T'.S = P.S xQ.S 
T'.i = (P.i, Q.i) 

T'.Tran = {((s P , s Q ),l P -l Q , (s' p , s' Q )) \ (s P , l P , s' p ) G P.Tran and (sq, Iq, s' Q ) G 
Q.Tran} 

T'.L = {I | (s,l,s') £ T'.Tran} 
T'.E = P.E U Q.E 
T'.I = P.I U Q.I 

We obtain T by discarding states in T' .S that are not reachable and discard- 
ing transitions in T'.Tran that are not satisfiable. We also remove any reference 
to events that appear in both P.E and Q.E — these are events that correspond 
to locations where two regions of excitable media meet, so they do not appear 
in the interface of P \\ Q. 

Definition 8. Forl P E P.L, Iq G Q.L, where I p -Iq is satisfiable, let 9(I p -Iq) be 
the expression obtained by setting any event variable a G P.EnQ.E that appears 
in l p ■ Iq to the value that leaves the resulting expression satisfiable. Constraint 
fj| and definition\Bp ensure that there will only be one such value for each a. 

Definition 9. The composition T = P || Q is defined as: 
T.S = {s G T' .S | s is reachable in T'} 
T.i = T'.i 

T.Tran — {(s,8(I p -Iq), s') \ (s,I p -Iq, s') G T'.Tran andl p -lQ is satisfiable} 
T.L = {k | (s, k, s') G T.Tran} 
T.E = T'.E \ (P.E n Q.E) 
T.I = T'.I \ (P.E n Q.E) 

We will now show that T = P || Q satisfies constraints [TJ [21 El and |3J 

Proposition 1. T satisfies constraint^ 

Proof. Take any ((s P , sq), A) G T.S x V(T.I). 

If P satisfies constraint 0] then let y p be the output subexpression corre- 
sponding to sp. Otherwise, by definition [6)3, Q.E n P.I = (i.e. no outputs 
from Q are inputs to P), and constraint [1] applied to P implies that there is 
exactly one (s p ,l p , s' P ) such that l p matches A n P.I, and let y p be the output 
subexpression of l p . 
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Similarly, if Q satisfies constraint 2] then let yg be the output subexpres- 
sion corresponding to sq. Otherwise, P.E n Q.I = and there is exactly one 
(sq, Iq, s'q) such that Iq matches AnQ.I, and let j/q be the output subexpression 
of l Q . 

Between them, the two subexpressions yp and uq determine a single C G 
V(P.EC\Q.E) for a given ((s P ,s Q ),A) G T.S x V{T.T). Constraint Q] applied 
to P and Q implies that there is exactly one (sp,lp, s' P ) such that lp matches 
A U C n P.I and exactly one (sq, Zq, s'q) such that Zq matches A\J C D Q.I. 

Therefore, for a given (s p , sq) there must be exactly one transition with label 
6(lp ■ Iq) that matches A. Thus constraint [T] holds for P \\ Q. 

Proposition 2. T satisfies constraint^ 

Proof. I is of the form v ■ w ■ x ■ y where v 6 Bool(P.I), w £ Conj(P.E \ P.I), x G 
Bool{Q.I) and y G Conj(Q.E \ Q.I). 

6(1) is obtained from I by setting all variables in P.EDQ.E to True or False, 
so 0(1) is of the form v' ■ w' ■ x' ■ y' where u' G Bool(P.I \ (P.E n Q.-E)), w' G 
Conj(P.E\(P.I U P.E HQ.E)),x' G Bool(Q.I \(P.E DQ.E)),y' G Conj(Q.E\ 
(Q.IUP.EHQ.E))). This can be written as (V -x^O' V)- Noww'-a:' G Bool(T.I) 
and w' • y' € Conj(T.E \ T.I), so constraint [5] is satisfied. 

Proposition 3. T satisfies constraint^ 

Proof. T was obtained it from T" by discarding those states that were not reach- 
able. 

Proposition 4. T satisfies constraint^ 

Proof. Since P and Q both satisfy constraint [21 for a given state (sp, sq) G T.S 1 , 
the label of every transition ((sp, sq),6(Ip ■ Iq), (s' p , s'q)) can be written as 
9(x P -yp-XQ ■ y Q ) = 9(x P ■ x Q ) ■ 9(y P ) ■ 9(y Q ). 

If P satisfies constraint [4] then yp is determined soley by sp (i.e. it is inde- 
pendent of s'p), and so is 9(yp). Otherwise, by definition [BJ3, all variables in yp 
are in P.E Pi Q.E, and so 9(yp) = True, so again 9(yp) is determined by sp. By 
the same argument, 9(yo) is determined by sq. 

Therefore 0(y p ) -9(yq) G Conj(T.E\T.I), is a subexpression common to all 
transitions beginning from (s p ,sq) G T.S, so T satisfies constraint [4j 

4 Examples 

4.1 Describing and composing fork and line structures 

We begin by considering single-use fork and line structures. Systems made from 
single-use structures behave like toppling domino systems. Each event can occur 
only once, and the system cannot return back to its initial state. In section 14.21 
we will deal with regions that recover after a period of time and which can be 
used again. 



Transition systems for excitable media 



17 



The transition systems given below for fork and line structures should be 
treated as axiomatic. No justification is offered here for why these transition 
systems describe the behaviour of fork and line structures: it is assumed that 
fork and line structures are sufficiently simple that the correspondence between 
their behaviour and the transition systems given below is obvious. However, once 
we have specified transition systems for the fork and the line, the behaviour of 
any structure built from these can be deduced automatically from the behaviour 
of the fork and the line by using the composition rules given above. 

Figure [TS] shows a unit length region of a horizontal channel in an excitable 
medium. The region has two input events a and c, and two output events b 
and d. Event a will cause event 6, event c will cause event d, but if both input 
events occur at the same time, there will be no output events. Figure \W\ shows 
a transition system that describes this behaviour. 




Fig. 15. A unit length horizontal channel in an excitable medium. 




Figure [T7l shows a unit length fork structure in the medium. Figure [18] shows a 
transition system that describes the behaviour of the fork, based on the informal 
description given in section^ 

To compose two of the excitable regions shown in Figure [TSJ let P be one 
region and let Q be a neighbouring region. Because we will be composing P and 
Q we would like event b of P to correspond to one of the input events of Q, and 
event c of P to correspond to one of the output events of Q, according to Figure 
[T51 We use the notation Q = P[m] to mean that Q is a transition system exactly 
like P, but with events renamed according to the mapping m. In this case we 
use the mapping m = {(a, 6), (d, c), (6, e), (c, /)} 
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Fig. 17. A unit length fork in an excitable medium. 




a,c,e:b,d,f 



Fig. 18. A transition system that describes the behaviour of the fork in Figure [TT1 




Fig. 19. Two composed unit length horizontal channels. 
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By definition [SI we can deduce P || Q from P and Q to obtain the transition 
system shown in Figure 1501 This transition system matches our intuitive under- 
standing of the behaviour of P \\ Q: event a will cause event e two time units 
later, unless event / occurs before the effect of a has time to propagate to e. 
Similarly, / will cause d unless it is interrupted by a. We can make a line of any 
length by composing unit length regions in this way. 



a 




Fig. 20. The transition system that describes the composition shown in Figure 1191 



Next we consider the single-use one way line shown in Figure |U Figure [5T] 
shows how this can be constructed from two forks F\ and F2, a line L\ of length 
5 and a line Li of length 3. L\ and Li are made by composing unit length 
regions, with their ends relabelled to match the events of F\ and Fi that they 
correspond to. 

Figure [55] shows the result of the composition F\ || L\ || Li || Fi- This 
figure introduces another notational convenience. Just as we use an arrow with 
an integer larger than 1 written above it to denote a scries of transitions of 
duration 1, so the split boxes in Figure [22l represent an interruptible sequence: 
a series of transitions that are followed so long as events that occur satisfy the 
expression inside the square brackets in the upper part of the box, but which 
can be interrupted on the occurrence of events that satisfy the expression in the 
lower part of the box. The number outside the square brackets in the upper part 
of the box is the length of the sequence. The arrow leading from the upper part 
of the box leads to the state that will be reached if the interrupting events do not 
occur, that leading from the lower part leads to the state that will be reached 
if they do. To illustrate this further, Figure 1531 depicts an interrupt box and the 
collection of states and transitions that it denotes. Sequences of states that fit 
patterns similar to this occur frequently in the type of medium considered here 
— for example an excitation propagating along a line of length n will reach the 
end after n units of time, unless interrupted by an excitation coming in the other 
direction. 
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Fig. 22. A transition system describing the behaviour of the one-way line Fi \\ L\ | 
L 2 || F 2 . 
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It can be seen from Figure [25] that the occurrence of a\ will lead to the occur- 
rence of &2 seven time units later (unless interrupted by C2), but the occurrence 
of C2 does not lead to any occurrence of d\ . 




a 



Fig. 23. The interrupt box notation above is equivalent to the system of states below. 



4.2 Reusable regions 

So far we have considered single-use regions. We can easily change a transition 
system describing a single-use unit length fork or line region into a transition 
system for a region that can be re-used after a 'recovery time' by adding a 
delay transitions from the final state back to the initial state. Figure |2"41 shows a 
transition system for the reusable equivalent of Figure (TS] with a recovery time 
of 5 time units. The choice of five time units is somewhat arbitrary, but leads 
to behaviour that is qualitatively similar to that obtained in BZ reaction media, 
where a propagating wave leaves behind it a short-lived inhibiting tail. The 
inhibiting tail prevents the propagation of further waves until it has disappeared. 

Let us consider a reusable one-way line, with the same structure as that in 
Figure 1211 but where each unit length region is reusable with a recovery time of 
5 time units. The transition system for this reusable one-way line has 41 states. 
The reason for such a large number of states is that as an event propagates 
through the system, different regions of the structure will recover at different 
times, and at any time after the inputs to the one-way line have recovered, 
further input events can occur. 

We can manage this extra complexity to some extent by placing constraints 
on the way in which a structure may be used. The context in which we are 
using a structure may provide us with knowledge about the number of times 
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a,c:b,d 



Fig. 24. A transition system for a unit length reusable region of an excitable medium. 



each input event can occur, or the rate at which they may occur. For example, 
we can specify that each input event must occur exactly once, in any order, 
with an arbitrarily long period of time separating the two events. This will not 
necessarily result in the same collection of behaviours as the single-use one way 
line, because we still allow the possibility that regions within the structure can 
be re-used. 

In order to specify this we can construct a second transition system R that 
restricts the pattern of input events to those that we are interested in, and then 
form T || R to produce a transition system limited to the behaviour determined 
byR. 

Let (Fi || L\ || L 2 || Fa)[{(ai, a[), (02, c' 2 )}} be the transition system for a 
reusable version of the one way line from Figure [2"T1 with a\ renamed to a[, and 
C2 renamed to c 2 . 

Figure [2"5l shows a transition system R with input events {oi, C2} and output 
events {a'i,c 2 } in which the first occurrence of a\ immediately leads to a' 1; and 
the first occurrence of C2 immediately leads to c 2 . Each of the output events a[ 
and c' 2 occurs exactly once in any path from state to state 2 in Figure HSJ R 
does not satisfy constraint 3J outputs a[ and c' 2 of R mirror the corresponding 
inputs immediately without any delay (for the first occurrence of each input). 
However, because R has no inputs connected to outputs of the system it is being 
composed with, and because all outputs of R feed into the system it is being 
composed with, the composition ((Fi \\ L\ || L 2 || F 2 )[{(ai, a[), (02, c 2 )}]) || R 
does satisfy constraint |U This composition is shown in Figure [Ml 

The reason why the rules for composability given in definition [5] are formu- 
lated to allow one of the transition systems being composed to flout constraint 
[?] is to permit us to do what we are doing here: they permit us to compose a 
transition system with a restriction transition system in which certain output 
events mirror certain input events. 

Figure [26] has several things in common with the single-use one-way line 
described by Figure [5U — in each case a\ will cause event b 2 if not interrupted 
by C2, but in Figure [55] the occurrence of event C2 before event b 2 is due to occur 
will always prevent 62 from occurring. In Figure l2"6l however, there is a window 
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^c2 




Fig. 25. A restriction transition system in which a[ and c' 2 both occur exactly once. 



al.c2 




Fig. 26. A transition system for the reusable one way line in which events a\ and C2 
occur exactly once. 
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of time during which the occurrence of ci after a\ will delay the occurrence of 62 
rather than prevent it. The transitions for this behaviour are shown as dashed 
arrows in Figure 1261 This is a behaviour of the reusable one-way line that the 
single-use version cannot exhibit. It arises from the possibility that regions within 
the structure can be excited more than once, even though the input events may 
not occur more than once. Figure [57] shows successive snapshots of a simulation 
of a sub-excitable BZ-reaction medium exhibiting this behaviour. 




Fig. 27. A one way line exhibiting the behaviour corresponding to the sequence of 
dashed arrows in Figure [26] 



4.3 A modulo-2 counter 

The final example considered here is a modulo-2 counter, shown in Figure |2"51 
This example is made from reusable excitable regions. It contains two one-way 
lines labelled 0\ and O2, seven forks labelled F\ to F7 and seven lines labelled 
L\ to L7. The length of each line is written in brackets in Figure [28] This 
example contains a loop structure made from L2, Fq, O2, L4, F4, F3, Lq and F% 
around which an excitation can propagate endlessly. The event c\ will lead to 
an excitation in this loop as follows: an excitation from C\ will split at fork F\, 
one excitation will propagate along L\ and into the loop travelling clockwise, the 
other will propagate through 0\ and into the loop travelling anti-clockwise. The 
latter will reach O2 and stop well before the former reaches O2 in the opposite 
direction, so the former will propagate endlessly around the loop. 

This is confirmed by the transition system in Figure 1291 which is produced 
from the parallel composition of all of the components, after specifying that no 
excitation will ever be input into the forks from which events /s and f? emerge, 
and that only one excitation will be fed into c\ . The fork F? provides an output 
fj which can be used to 'see' the excitation every time it propagates around the 
loop. 

Every time around the loop, the excitation will split at Fq and one excitation 
will travel along L$. It will not reach fa because it will be annihilated before it 
reaches F§ by an excitation splitting from the loop at P4. 

Once there is an excitation in the loop, it will periodically progress along L3 
into 0\. It will be stopped by 0\. If a second input is applied to ci at the right 
moment, it will coincide with the hrst excitation in the loop, so will not cause a 
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f q 

Fig. 28. A structure that can behave like a modulo-2 counter. 



-.cl 




cl:bl,f5,f7 



Fig. 29. A transition system for Figure [55] when only one input excitation is applied. 
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second clockwise-travelling excitation. It will also meet the first excitation in 0± . 
Its progress through 0% will effectively be delayed, according to the behaviour 
of 0\ given by the dashed line in Figure [26l so it will arrive at O2 late enough 
to cause the first excitation to be delayed in its progress around the loop, also 
according to the behaviour in Figure [211 This will allow the excitation travelling 
along L5 to reach output /s, and also to reach the loop and annihilate the 
excitation in the loop. 

To summarise: a single excitation at c\ will cause an excitation to propagate 
around the loop. A second excitation at c\ (at the right moment) will lead to 
an output at f§, and will return the system back to its original state. Therefore 
an output will emerge from /s for every other event at c\ (subject to timing 
constraints). This is the behaviour of a modulo-2 counter. 

Inspection of the 191-state transition system (not shown here) for the case 
when ci occurs twice reveals that in order for the system to exhibit this be- 
haviour, the interval between the first and second occurrences of c\ must be 
35,36 or 37 time units. By connecting a transition system that produces an out- 
put every 36 time units to input c\ , the transition system shown in Figure 1301 is 
obtained. This shows that, as expected, there is one occurrence of f$ for every 
two occurrences of /V, and that the system reverts to a previous state after the 
occurrence of f§. 




:bl,f5,f7 



Fig. 30. A transition system for Figure |2"51 when an input excitation is applied every 
36 time units. 



5 Implementation 

An implementation of a set of routines for manipulating transition systems, 
written using the Python 3 programming language, is available from: 



http : //www. srm. org.uk/downloads/transition.py 



This source file contains functions for creating all of the examples that are 
used in this paper. The Python programming language was used because it has 
built in data types for representing sets, tuples and maps, because it has a concise 



Transition systems for excitable media 



27 



syntax, and because it is freely available and widely supported. The source file 
contains comments which will help with understanding the implementation. 

6 Discussion and conclusions 

This paper has demonstrated that it is possible to use transition systems to 
describe the behaviour of unit length fork and line structures in excitable media, 
and that it is possible to deduce the behaviour of larger structures from the 
behaviour of their components. Although this paper has exclusively used fork 
and line structures for making systems, it should be possible to use the method 
given here for automatically deducing the behaviour of other excitable media 
structures made from simple components, where the behaviour of the simple 
components can be described by transition systems. 

I have chosen to model behaviour that involves the propagation and inter- 
action of excitations in channelled structures, but there are other types of be- 
haviour that excitable media can exhibit, and which can be used for information 
processing, that the methods used in this paper cannot describe. The use of a 
planar BZ-reaction medium for image processing [25] exploits a large rectangular 
section of the medium and so cannot be described in terms of excitations propa- 
gating along channels. Steinbock, Kettunen and Showalter [26] make use of the 
interaction of several waves over an extended area to implement logic gates. A 
transition system could be used to describe one of these logic gates and its tim- 
ing properties, but it could not be used to describe the lower-level extended-area 
interactions of excitations and deduce the behaviour of a logic gate from these. 

Much research has been carried out on the use of process algebra to describe 
and reason about concurrent systems. There are two main reasons why I chose to 
work directly with transition systems and use a graph-based representation for a 
set of behaviours, rather than make use of process algebra and adopt a linguistic 
representation for a set of behaviours. Firstly, process algebraic notation can 
often obscure the simplicity that is readily apparent in a graphical representation: 
even for small systems, not all sets of behaviours can be represented in an easy- 
to-understand way in any of the current linguistic representations employed in 
process algebra. Secondly, the range of meanings and structures that can be given 
to labels in transition systems provides enough flexibility to permit the physical 
systems considered here to be described, but the structure of transition systems 
is constrained enough to reduce the danger of losing the ability to compose 
systems together. 

However, we have only considered relatively simple systems in this paper and 
I do not claim to have fully explored all of the possible ways of using concurrency 
theory to model structured excitable media. For more complex systems a well- 
formulated algebraic approach that permits equational reasoning and recursive 
definitions may be useful. The sequence and interrupt box notations in Figures 
[TTl and [221 respectively were introduced as syntactic conveniences for represent- 
ing collections of states and transitions. Instead of regarding these simply as 
syntactic entities they could alternatively be treated as semantic entities, and 
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rules for manipulating them during composition could perhaps be formulated, 
without having to expand them into collections of transitions and states. 

Several simplifying assumptions and approximations have been made. I have 
assumed that all regions of the medium that we are modelling propagate excita- 
tions at the same rate. This assumption is valid for simulated planar BZ-reaction 
media. It is also reasonable in experimental planar BZ-reaction media, where ex- 
perimental conditions can easily be achieved in which the propagation of a wave 
of excitation through two identically structured regions takes a very similar 
length of time. Theoretical and experimental investigations of domino lines indi- 
cate that their behaviour is also predictable and repeatable |27] . This assumption 
is less valid for plasmodium of physarum polycephalum, where experiments to 
measure the speed of propagation of plasmodia in the spreading/contracting 
mode [22] show that the speed of propagation can vary over the course of a sin- 
gle experiment, and can vary even more from one experimental run to another. 
It is not yet known how easily these variations can be reduced by controlling the 
experimental conditions. 

I have used discrete time rather than continuous time. The effect of this 
is to restrict lengths in structures to integer multiples of the length that an 
excitation covers in a single time unit, and to restrict the times at which input 
events can occur to integer multiples of a single time unit. Since we are free to 
choose which units of physical time and length to use, this does not seem to be a 
severe restriction. It is unlikely that it would be possible to make a deterministic 
physical structure so critically dependent on timing for its operation that it 
cannot be modelled using integral channel lengths for some choice of length unit 
and time unit. 

Modelling a physical system using a transition system involves identifying the 
states of the physical system with states of the transition system, and identifying 
changes of state in the physical system with transitions in the transition system. 
Trajectories of a system correspond to sequences of states and transitions. For 
synchronous systems where each trajectory has a common clock, an alternative 
approach can be conceived of in which trajectories rather than states are treated 
as primary, and transitions in a system are from one trajectory to another rather 
than from one state to another. If such a scheme can be formalised, and if parallel 
composition rules for such a trajectory-oriented scheme can be formulated, then 
this might be a more natural way of modelling the type of system considered in 
this paper. 
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